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ABSTRACT 

The first macroscopic bodies in protoplanetary disks are dust aggregates. We report on a number of experimental studies with dust 
aggregates formed from micron-size quartz grains. We confirm in laboratory collision experiments an earlier finding that produc- 
ing macroscopic bodies by the random impact of sub-mm aggregates results in a well-defined upper-filling factor of 0.31 ± 0.01. 
Compared to earlier experiments, we increase the projectile mass by about a factor of 100. The collision experiments also show that 
a highly porous dust-aggregate can retain its highly porous core if collisions get more energetic and a denser shell forms on top of 
the porous core. We measure the mechanical properties of cm-sized dust samples of different filling factors between 0.34 and 0.50. 
The tensile strength measured by a Brazilian test, varies between 1 kPa and 6 kPa. The sound speed is determined by a runtime mea- 
surement to range between 80 m/s and 140 m/s while Young's modulus is derived from the sound speed and varies between 7MPa 
and 25 MPa. The samples were also subjected to quasi-static omni- and uni-directional compression todetermine their compression 
strengths and flow functions. Applied to planet formation, our experiments provide basic data for future simulations, explain the spe- 
cific coUisional outcomes observed in earlier experiments, and in general support a scenario where coUisional growth of planetesimals 
is possible. 
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1. Introduction 

Collisions have an important influence on the evolution of plane- 
tary systems from their early beginnings onward. Understanding 
the effect of a collision requires knowledge about the mechanical 
properties of the colliding objects of given sizes. The different 
size ranges considered for the very early phases of planet forma- 
tion are frequently termed pebbles, rocks, or boulders, though 
the mechanical properties of the initial population do not reflect 
these names. Initially, all objects can be considered as dust ag- 
gregates. As bodies smaller than planetesimal size (~ km-size) 
couple sufficiently well to the gas of a protoplanetary disk, col- 
lision velocities stay moderate and 'only' reach several tens of 
m/s ( Weidenschilling & Cuzzi|[l993| l. At these velocities, col- 
lisions are in general insufficiently energetic to melt or petrify 
an object. Therefore, originating from dust particles in a proto- 
planetary disk, the initial population of objects have to be dust 
aggregates of a wide range of sizes. 

Two different approaches currently exist to simulate colli- 
sions between dust aggregates. The first method is based on 
molecular dynamics and simulates the interaction between each 
dust grain in two colliding aggregates. This method is restricted 
to small (<; 1 mm) dust aggregates owing to computational 
limits (|Dominik & TielenslfTW?! |Wada et al.|[2009l [Seizinger 
|et al.|[ 2012[ l. The other approach uses smoothed particle hy 
drodynamics (SPH) to simulate dust aggregates (Schaf er et al.| 
[2007, ,Benz"& Asphaug||1999[ [Sirono„200 4). The single grains 



are no longer resolved in this method, but aggregates are charac- 
terized in terms of macroscopic mechanical properties, such as 
their compressive strength, tensile strength, or Young's modulus 
( [Geretshauser et al.|201 la| l. The SPH-method is not restricted to 
a certain particle size, hence can provide crucial input for global 



coagulation models. However, the mechanical properties of pro- 
toplanetesimals are first needed as input for the simulations. One 
way to obtain these parameters is to measure them in laboratory 
experiments. 

So far, only a few studies of the mechanical properties of dust 
aggregates have been performed. Blum & Schrapler (200 4) per- 
formed the first experiments to produce and characterize highly 
porous dust-aggregates (produced by random ballistic deposi- 
tion) with a filling factor of 0.15 for a dust material consisting of 
monodisperse, spherical silica grains (d = LSj-im). They deter- 
mined the evolution of the filling factor during an uni-directional 
compaction of the dust aggregates. In a later study, these ex- 



periments were extended to different dust materials dBlum et al. 



2006). In both of these studies three different dust samples were 
used (monodisperse silica spheres, diamond powder, and poly- 
disperse quartz) and the porosity evolution at increasing uni- 
axial pressure was determined. Additionally, the tensile strength 
was measured for highly porous samples, produced by random 
ballistic deposition. [Giittler et al. ( 2009| l determined the omni- 
directional compression, but only for the monodisperse dust 
sample previously studied by Blum & Schrapler (20041 and 
Blum et al. ( |2006 i. The mechanical properties strongly depend 
on both the shape and size of the monomers. This was shown ex- 
perimentally by Blum et aL] ( |2006) and theoretically by [Bertini 
|etaL] ( |2009l l. 

The experiments presented here are dedicated to providing 
the mechanical properties for a polydisperse dust material with a 
grain size distribution, which might be typical of protoplanetary 
disks. Different experimental methods are applied to measure 
critical parameters such as tensile strength, sound velocity, or 
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compressive strength. The mechanical properties are studied for 
various filling factors. 

Besides determining the basic properties of small dusty bod- 
ies, this work contains experimental results on porosities that 
might be expected for small sub-meter size protoplanetesimals. 
Further experiments also show that an evolving protoplanetesi- 
mal is unlikely to be homogeneous but we show that core - shell 
particles with highly porous cores and more compact shells can 
easily form through collisions with small mm-size aggregates. 

The paper is organized as follows. We first briefly describe 
the dust samples and their preparation as generally as possible 
in section 2. We then give details of the experiments carried out 
with these samples in section 3. Results and some general in- 
terpretation are given in section 4 and the final section 5 places 
the results and conclusions in the context of application in pro- 
toplanetary disks. 



2. Sample preparation 

For all experiments, we used quartz dust with particle sizes rang- 
ing between 0.1 /im and 10 jum (80% of the mass consisting of 
particles of size between 1 fj.m and 5 yum) from which we pro- 
duced larger cm-sized aggregates. The same quartz dust had 



been used in several previous experiments (Blum et al. 2006 
IWurm et al.|2005al[Tiiser & Wurm|2009l l. We shidied a number 
of different properties of the dust aggregates depending on the 
(volume) filling factor 



V. 



solid 



V, 



(1) 



total 



which is the ratio of the volume filled by solid material (VsoUci) 
to the total volume (V/om/)- This is equivalent to the ratio of 
the density of the total object to the density of its solid quartz 
building-stones, which is 2.6g/cm''. The samples for further 
measurements were always formed by compression of dust pow- 
der from all sides as the sample was confined to a cylindrical 
mold. Pressure was applied to the top face of the cylinder, which 
was free to move (Fig.[T^). During the preparations, the applied 
omnidirectional pressure cri was measured using a force sen- 
sor attached to the compressing piston (Fig. [T^). After correc- 
tion for the exerted pressure to the additional weight of the top 
face bolt, ctj represents the compression stress that the dust ag- 
gregate was exposed to (Fig. [T]?). The maximum pressure ap- 
plied was 55 kPa. We prepared cylinders with constant masses 
of 20 g that always had a diameter of 30 mm. For this reason, 
their heights varied between 22 mm and 36 mm for different fill- 
ing factors. After removing the dust cylinders from the mold, the 
volumes and detailed masses were measured. With this informa- 
tion, densities and therefore filling factors could be determined. 
In total, the filling factors of the analyzed dust samples ranged 
from 0.24 to 0.51. The presented values for filling factors are 
mean values for whole samples. For the more porous samples 
in particular it is unclear a priori whether the dust cylinders are 
homogeneous. To estimate any inhomogeneities, we produced 
dust cylinders with different filling factors in the range and way 
described above. We cut them into two halves with respect to 
height and measured their filling factors in the same way. Each 
top half, which was more affected by the pressure applied to the 
top face of the cylinder was slightly more dense than the bot- 
tom half The difference between the bottom and top halves was 
about 0.02 in filling factor, independently of the calculated mean 
value within our investigated range. These inhomogeneities - 
which are ±0.01 deviations from the mean - are small compared 



to the overall variations in the filling factor and are neglected in 
the following analysis. 



3. Experimental setups 

The following different experiments were carried out to specify 
the elastic properties of the dust aggregates and some general 
trends in the evolution of the filling factor. As far as the mechan- 
ical properties are concerned, we studied: 

- omni- and unidirectional compression, 

- tensile strength, 

- sound velocity, 

- modulus of elasticity (Young's modulus) based on the sound 
velocity measurements. 

Besides the determination of these basic properties, collision ex- 
periments were carried out to show the influence of collisions on 
the evolution of the filUng factors of a forming dust aggregate. 



3. 1 . Procedure for compression measurements 




Fig. 1. Experiment to measure the compression strength of a 
dust sample: a) production of a cylindrical dust- aggregate by 
applying a certain solidification pressure to a confined dust sam- 
ple; b) removal of dust cylinder produced by a given compres- 
sion stress (Ti; c) the compression strength cr^ is measured when 
the sample breaks ( |Jenike|1964| l. 

Cylindrical dust-aggregates were produced as described in 
section |2] To these cylinders, we applied an unidirectional pres- 
sure at the top face and determined the threshold pressure cTc 
(equivalent to compression strength) when the sample did break 
as sketched in Fig. [T];. We corrected the exerted pressure to the 
additional weight of the piston and obtained the uniaxial com- 
pression strength for each precompacted dust cylinder. Samples 
were used with different precompactions i. e. filling factors rang- 
ing from 0.24 to 0.42. 

3.2. Procedure for tensile strength measurements 

To determine the tensile strength of the dust cylinders, we per- 
formed Brazilian tests. We note that independent but apparently 
similar measurements have been carried out by Kothe et al. (per- 
sonal communication) for the same type of dust and it will be 
interesting to perform a detailed comparison with their results in 
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the near future. The essence of the BraziHan test is that a pressure 
is imposed upon the lateral surface of the dust cylinders (Fig.|2]i 
in terms of two opposing effective single-loads. The pressure at 
which the sample breaks is then determined. 




Fig. 2. Setup for the Brazilian test: a piston connected with a 
force sensor is pushed down onto a bar lying above the sample. 
The two bars are used for a linear load distribution. The brackets 
hold the cylinder in place initially. Pressure is applied until the 
aggregate breaks. 



A linear distribution of an applied load must be provided as 
illustrated in Fig.[3j to prevent too large compression stresses at 
the point of load application. It is necessary to produce a tensile 
fracture when the dust cylinder is crushed. According to Mitchell] 
( 1961| l, a fracture of a breaking dust cylinder in a splitting ten- 
sile test can take different shapes depending on the width of the 
load distribution bars. While the non-use of bars causes a stress 
fracture, only the use of narrow bars generates a vertical division 
line characteristic of a proper fracture, which is the indication of 
a tensile fracture. In practice, two small wedged pieces of dust 
are produced at the break. They develop directly above and be- 
low the bars and indicate a shear fracture. We used bars with 
a width of 5 mm and a height of 2 mm. Single vertical fracture 
lines were always produced from nearly the top bar to the bottom 
bar. The detailed width or contact area does not enter into the fi- 
nal determination of the tensile strength as only two times the 
applied force and the mantle surface of the dust cylinder enter 
according to 



for diameter d and length L of the cylinder. 



load 




Fig. 3. Illustration of the load distribution in a dust cylinder dur- 
ing the performance of the Brazilian test. The fracture pattern is 
given along the dashed lines according to |Mitche"il| ( | 1 96 1 | l . 



The bars were placed between the dust cylinder used and the 
acting loads from the piston and the table. The dust aggregate 
was confined between two brackets to avoid rolling. On top, a 
piston and a force sensor were attached to one bar to determine 
the splitting force. The two brackets were removed as soon as 
the piston is in place and the load sufficient to keep the dust 
cylinder in place. Pressure was then applied. In addition to the 
pressure measurements, we used a high speed camera to observe 
the formation of the fracture and the break-up of the samples. 
For these experimental series, samples were used for a range of 
volume filling-factors from 0.37 to 0.51. 

3.3. Procedure for sound velocity measurements 

The sound velocity in the dust cylinders was measured in a basic 
way. A sound wave was initiated through a mechanical contact 
on one side of the aggregate by manually knocking against the 
aggregate support structure. Our results do not depend on any 
variations in the manual knocking. As the sound velocity was 
known to be lower than the speed of sound in air, we measured 
the sound velocity in a vacuum chamber below 4- 10"^ mbar (Fig. 
|4|l. The time when the pulse is generated was recorded, the pulse 
itself closing an electrical contact. The arrival time of the sig- 
nal at a force sensor on the other side of the aggregate was also 
recorded. The time delay was corrected for sound-wave travel 
times using the attached hardware. With the predetermined ag- 
gregate size, the sound velocity was calculated. We averaged 30 
to 40 measurements for a given aggregate. These measurements 
were performed with dust cylinders that varied in terms of their 
filling factors across the range from 0.33 to 0.51. 

3.4. Procedure for collision experiments with dust particles 

The previous measurements were restricted to determining the 
properties of dust aggregates that were artificially formed with 
a defined volume filling-factor resulting from the omnidirec- 
tional compression within the mold. With protoplanetary growth 
in mind, compression might however result from collisions. A 
likely impact in the early formation phases is the impact of a 
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vacuum 
chamber 



beat generates 
sound wave 



dust 

aggregate 




force sensor 



Fig. 4. Setup for measurement of sound velocities in dust sam- 
ples inside a vacuum chamber 



smaller body (projectile) with a larger body (target). This impact 
can be described by neither uni- nor omnidirectional static com- 
pression. Building a dust aggregate from scratch by impacts of 
sub-mm dust aggregates, [Teiser et al. (201 la i found that homo- 
geneous aggregates form with a volume filling-factor that only 
depends on the impact velocity and reaches a limiting filling- 
factor of 0.31 at an impact velocity of 6 m/s. We used the same 
setup here. Targets were placed in a vacuum chamber at the 



and produced dust particles with a mean size of 250 fim. By us- 
ing a sieve with a mesh size of 4 mm, we were able to produce 
dust projectiles with a mean particle size of 1.2 mm. This is an 
increase in the mass of the projectiles by about a factor of 100. 

In a second set of experiments, we first produced highly 
porous targets through impacts at low velocities. We formed 
these highly porous targets by ejecting dust through the sieve 
of mesh size 125yL(m directly held above a transparent box at 
normal air pressure. We then raised the sieve to increase the im- 
pact velocity on the prepared highly porous target. In different 
experiment runs, we varied the sizes of the impacting fast dust- 
aggregates by using sieves with mesh sizes of 125 /im and 4 mm. 
We marked the boundary of the porous target with a colored 
layer of dust to illustrate the depth to which impacts influence 
the target below. 

4. Results and interpretations 

In the following we present our results and provide some inter- 
pretations. 



4. 1 . Compression measurements 

In the case of compression during sample preparation. Fig. 
|6] (top) shows the density and volume-filling factor of the dust 
cylinders over compression stress applied to the confining mold 
for production of the cylinders. By fitting a power law 

p=Po-(l + — )" (3) 

characterizing the isentropic compression of bulk material, we 
can classify the SiOa dust used by a calculated index of com- 
pressibility n ( [TomaspOOOj |2001 !). Additionally, we obtain the 
isostatic tensile strength (Tq and the initial density po of the dust 
placed in the mold. The following values were determined: the 



Table 1. Parameters for isentropic compression of quartz dust 
samples 



transparent 
box 




dust particles 
V = up to 7m/s 



very porous 
dust target 



Fig. 5. Dust aggregates are generated by impacting dust projec- 
tiles with velocities up to 7 ms"'. 



bottom of a tube. After evacuating the whole tube, dust parti- 
cles of different sizes were ejected through a sieve placed in the 
tube at different heights. Accelerated by gravity in free fall, dust 
aggregates impact the target with predefined velocities between 



1 ms and 7 ms (F: 



; targ( 



In a first set of experiments, we ex- 



tended the study of Teiser et al. (201 la i to larger projectiles. 
|Teiser et aLlpOlla] ) used a sieve with a mesh size of 125/im 



Po (g/cm^) 


o-o (Pa) 


n 


0.612 ±0.024 


295.26 ± 155.70 


0.115 ±0.009 



index of compressibility characterizes the compressibility and 
cohesion / flow ability according to the distinction of table [2] 



( TomaspOOl i. In our case, the dust cylinders are classified as 



very compressible and very cohesive. 

Table 2. Semi-empirical classification of the index of compress- 
ibility 



index n 


evaluation 


flow ability 


0-0.01 
0.01 - 0.05 
0.05-0.1 
0.1 - 1 


incompressible 
low compressibility 

compressible 
very compressible 


free flowing 
(free) flowing 

cohesive 
very cohesive 



In the case of unidirectional compression, an additional 
classification can be based on the flow function ffc, which is de- 
fined as the ratio of the compression stress cri to the compression 
strength of the samples cTc 



0"c 



(4) 
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Fig. 6. Top: density and filling-factor dependence of quartz 
dust on the compressive stress; Bottom: uniaxial compression 
strength as a function of the omnidirectional compression stress. 
Values for the flow function cr^ - i{cr\) can be extracted from 
the intersections of the dashed lines with the data (red line). 

The flow function is illustrated in Fig. [6] (bottom), which shows 
the compression strength versus compression stress. The data 
show a typical trend, as more solidified dust cylinders have 
greater resistance to the applied unidirectional pressure. The 
flow function leads to the following classification as given in 
table [3] ( Jenikel 1 96 4*). The characterization varies between "not 
flowing" for ffc < 1 and "free flowing" for ffc > 10: We find that 



Table 3. Flow ability of granular matter or powders (Jenike 
[T964l l 



fF, < 1 


non flowing 


1 < ffc < 2 


very cohesive 


2 < ff, < 4 


cohesive 


4 < ff, < 10 


easy flowing 


10 < ff. 


free flowing 



low density cylinders (low applied compression stress) fall into 
a class of being cohesive or very cohesive in agreement with the 
classification given above. In contrast to that, the most compact 
dust-aggregates are characterized as flowing. This does not con- 
tradict the classification from the power-law fit for the isentropic 
compression given above, as this is a single fit to the whole range 
of compression stresses but is mostly sensitive to the low com- 
pression stress values where the density changes significantly. If 
the behavior of the dust aggregates were also divided between 



low and high compression stress, the data for the high compres- 
sion stress would be in agreement with a much lower index of 
compressibility and therefore flowing conditions (Fig. |6] top). 
The flow function seems to be a more reliable tool here to char- 
acterize dust aggregates of different filling-factors. 

The resulting classification then also agrees with general ob- 
servations that aggregates of either low density or low filling fac- 
tor can more easily lose energy in collisions by restructuring the 
constituent grains and compaction when dense aggregates are 
more elastic ( |Beitz et al.|2011p . 

4.2. Tensile strength measurements 

During the Brazilian test shown in Fig. |2j dust cylinders with 
different densities were observed to crack into two halves. An 
example image sequence is shown in Fig. |7] Additionally, two 




Fig. 7. Image sequence of the disruption of a dust cylinder into 
two halves during a Brazilian test. The expansion of the fracture 
from the middle towards the bars is clearly visible. 



very small wedged pieces of dust caused by the use of the load 
spreading bars, were produced at the break. They developed di- 
rectly above and below the bars. According to |MitcheIl| ( | 1 96 1 | l, 
this is an indication of a shear fracture. We also observed that 
the fracture developed in the middle of the cross-section and 
propagated in the direction of the two bars. This is in agreement 
with the investigations of ^Rocco et al. ( 19991. The dependence 



of the splitting tensile strength on the volume filling-factor is 
seen in Fig. [8] With an increase in volume filling-factor <p, the 
dust aggregates show more resistance or have larger values of 
the splitting tensile strength crtj. Without seeking further justifi- 
cation but describing the increase in the tensile strength in terms 
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40 



0.3 



0.35 



0,4 0.45 
volume filling factor (|) 



0.5 



0.55 



Fig. 8. Splitting tensile strength over density measured for 21 
dust cylinders (open circles). The data were binned (black 
crosses). The red line shows an exponential fit. 

of increasing filling-factor, the most convenient fit turned out to 
be the exponential dependence 

,14.3-0 



cr,s = 4 ■ e^ 



(5) 



4.3. Sound velocity measurements and modulus of elasticity 
(Young's modulus) 

In Fig. [9] we present our measured data for the sound velocity. 
With increasing filling-factors, a linear increase in the sound ve- 
locity Vs can be detected, which we fitted as 

Vs = 337 ■ 0m/s - 26.4 m/s. (6) 

The absolute values are small in comparison with sound veloc- 
ities in monolithic quartz (vj x 5000 m/s). The obtained sound 
velocities in powdery Si02 samples are several 10 times smaller. 
We note that the comparison to measurements not done in a vac- 
uum but at an ambient pressure of 1 atmosphere showed no sig- 
nificant difference. 

As the modulus of elasticity is not measured easily for 
porous dust-aggregates, the sound velocity provides a means for 
an indirect determination. The sound velocity is 



Cl 



E(l-y) 



|p(l+v)(l-2v)' 



(7) 



where E is the modulus of elasticity, p is the density, and v is the 
lateral contraction or Poisson number. By neglecting the Poisson 
number, which is usually small (v ^ 0), and solving the equation 
for E, we obtained the estimation for our modulus of elasticity 

E^Ci^-p (8) 

To calculate the moduli of elasticity, we performed a linear 
fit for the sound velocity and identified Vj with Ci. Together with 
p = 2600 kg/m^ ■ 0, we derive 



E = 2.96 ■ 10*^ ■ (p^ - 4.68 ■ 10^ ■ 0^ + 1.82 ■ 10^ ■ </> ■ 



N 



m^ 



(9) 



10 The 



The resulting moduli of elasticity are shown in Fig. 
obtained data are very low values for the moduli of elasticity 
or Young's modulus (on the order of 10 MPa). The modulus of 
elasticity for soUd SiOa is given by 85 GPa ( jCarlotti|1995| l which 
is about three orders of magnitude larger. 



0.3 



0.35 



0.4 0.45 
volume filling factor i(l 



0.5 



0.55 



Fig. 9. Measured data for sound velocity versus filling factors. 
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Fig. 10. Calculated moduli of elasticity based on the measured 
sound velocities. 



4.4. Results and observations of collision experiments with 
dust particles 

4.4.1. Homogeneous aggregate formation 



[Teiser et al. ( |2011a i used particle sizes of 0.25 mm to grow dust 
aggregates. In this study, we used a larger mesh size for the sieve. 
This provides dust particles with mean particle sizes of 1.2 mm. 
Approximating the mass by the third power of the size, this is a 
factor of 100 more in mass. The dependence of the filling-factor 



on the impact velocity is compared in Fig. 1 1 with the data mea- 
sured by Teiser et al. ( 201 la I. At lower velocities (approximately 
up to 3 m/s), the generated dust aggregates are more compact if 
built from larger particles. This is understandable as they pos- 
sess more energy and can restructure the target surface more ef- 
ficiently. However, once the collision velocity reaches 5 m/s, the 
forming aggregates reach the same limiting volume filling-factor 
of 0.31 ±0.01. 



4.4.2. Layered aggregates 

As basic targets in this set of experiments, we used highly porous 
dust targets formed by slow impacts of aggregates produced of 
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0,15 



mean particle sfze 0.25mm, Teiser et al. 2011 
mean particle size 1.2mm 



■ f- 



3 4 5 

impact velocity [m/s] 



Fig. 11. Volume filling-factors of dust aggregates generated by 
free falling impacting dust aggregates for two different projectile 
sizes; blue 1.2 mm, this study; red 0.25 mm, from |Teiser et al. 
( |2011a| l 



0.25 mm aggregates (Fig. 12 top). These highly porous targets 
still show the characteristic size of the projectiles in their tex- 
ture. The initial aggregates (Fig. 12 la and 2a) have a mean 



volume filling-factor of 0. 14. We then exposed these targets to 
impacts at higher velocities of about 6 m/s. As aggregate size, 
we used in the first case (Fig. 



12 



lb) the small 0.25 mm and 
in the second case the large 1.2 mm particles (Fig. 12 2b). The 
small projectiles build up a compact layer above the initial target 
sample. The initial target is unchanged below the marked layer 
Only the initial, about 1 mm thick, marked layer of the target is 
therefore influenced by the impacts (Fig. 12 lb). A more ex- 



tended compression zone of about 6.5 mm was observed for the 
larger projectiles (Fig. 12 2b), even thoughthe initial target re- 
mained unchanged below. A quantification of the transition layer 
depending on velocity is currently beyond the capabilities of the 
setup. However, the experiments show that a highly porous dust- 
aggregate of cm size or larger can keep this porosity in its core if 
collision velocities increase and a more compact shell is added. 
Core mantle aggregates would be the result with inhomogeneous 
filUng-factors. 











5nnnn 













5. Application and conclusions 

The collisional evolution from //m-sized dust grains to km-sized 
planetesimals and beyond strongly depends on the mechanical 
properties of dust aggregates. The basic mechanical properties 
are the modulus of elasticity, the tensile strength, and the com- 
pression strength. These parameters are intrinsically difficult to 
determine for dust aggregates. It is difficult to measure the mod- 
ulus of elasticity directly, as a small amount of stress compresses 
the aggregates, changing both the porosity and modulus of elas- 
ticity. However, by measuring the speed of sound, which is re- 
lated to the Young's modulus, we were able to determine the 
porosity dependence of a dust sample. We also succeeded in ap- 
plying the Brazilian test to dust aggregates of moderate density 
and determining the tensile strength of the same kind of dust ag- 
gregate. Determining the compression strength in addition gives 
a set of properties, which allows us to develop a more realis- 
tic model for dust aggregates in numerical simulations of colli- 
sions, though the set is still incomplete as e.g. shear is not con- 



Fig. 12. Dependence of the volume filling factor on impacting 
projectiles; la, 2a; low filling factor aggregates (p - 0.14 gen- 
erated by sub m/s impacts; lb; change of original target upon 
collisions with 0.25 mm aggregates at 6 m/s; 2b; change of orig- 
inal target upon collisions with 1.2 mm aggregates; The transi- 
tion zone is only on the order of up to several mm. It separates 
two protoplanetesimal zones with filling factors characterized by 
the process by which they formed, i.e. a core-shell particle with 
highly porous core and dense shell forms. 



sidered. In total, different conclusions can be drawn from our 
experiments, which are listed first and detailed afterwards. 

- Most basicly, the results can be used as input parameters 
to detailed numerical simulations of collisions. In this re- 
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spect, the measurements allow no further conclusion here. 
Implications have to await their use in such simulations. 

- The porosity evolution in a self-consistent growth scenario 
can be compared to earlier speculations by collecting small 
dust aggregates at a given velocity. 

- The porosity is important as it strongly influences the ther- 
mal properties of dust aggregates and our results can be used 
to quantify possible consequences, e.g. for photophoretic 
particle sorting in protoplanetary disks. 

- Even without detailed simulations, the outcome of published 
collision experiments between cm-size aggregates can be re- 
interpreted in the light of the determined properties, depend- 
ing on possible porosity evolutions. 

- Sound propagation in cohesive and granular matter are part 
of the research in communities apart from planet formation 
and our measurements can be compared to more fundamen- 
tal models. The exact value is also important in helping us to 
decide whether a collision is supersonic or subsonic, which 
has an impact on the way in which we model dust-aggregate 
collisions. 

As noted, not much more can be said about the first item, which 
is an input parameter to numerical simulations. It is known 
that variations in material properties can strongly influence the 
outcome of a collision. Depending on either the porosity or 
the existence of a less porous (hard) shell on a porous core, 
two aggregates might stick, bounce off each other, or fragment 
( |Geretshauser et al.|20I lb||20I0||201 l a)). The full importance of 
our measured data in this range of outcomes can only be seen 
in detailed simulations that remain to be carried out in the fu- 
ture. Typical values of the parameters studied are in the follow- 
ing range: 

- compression strength (unidirectional) 1-8 kPa 

(for omnidirectional compression and flow function see sec- 
tionlTT), 

- tensile strength 1-5 kPa, 

- modulus of elasticity 5-25 MPa, 

- sound speed 80 - 140 m/s. 

Nevertheless, the measured properties already shed some 
light on recent collision experiments. Beitz et al. ( 201 l| l stud- 
ied the collisions between cm-size dust aggregates, which had 
porosities in the same range as the porosities studied here and 
the same kind of dust (micron-sized quartz). They found that in 
a collision at about 1 ms"', a slightly more porous aggregate gets 
destroyed and adds some mass to the more compact aggregate, 
which stays intact even if the porosity difference is only 0.02. 
This can be qualitatively explained by having a look at the ten- 
sile strength (Fig. [8]l, which shows an exponential dependence 
on the volume filling-factor. Therefore, a slightly more compact 
aggregate is much more stable. It was not expected that this ef- 
fect be so strong nor that very small variations in porosity could 
decide whether an aggregate fragments or growth. This clearly 
shows that the porosity evolution during planetesimal formation 
is a key factor in the process of planetesimal formation. 

While they were meant to be indirect tools for determining 
the modulus of elasticity, the speed of sound measurements are 
also important in their own respect. They vary between 80 ms"' 
and 140 ms ' across the range of studied porosities. This is 
much lower than the speed of sound in monolithic solid material. 
Numerical simulations by Paszun & Dominik ( 2008 ) and Ringl 



grains. For moderate models of protoplanetary disks, the abso 
lute values do not differ considerably but are clearly ab ove the 



collision velocities of up to 50 ms ' (Weidenschilling & Cuzzi 



1993). However, somewhat more turbulent disk collisions might 
become supersonic rather early, which can influence any subse- 
quent collisions. 

A few more aspects of collisional evolution come from our 
collision experiments in this work. They show that the average 
properties of a whole aggregate are appropriate if this aggregate 
grows by a continuous process, i.e. by adding small aggregates at 
a given collision velocity. The volume filling factor depends on 
the collision-velocity but we note that has a limit of about 0.31 
at high collision velocities. We extended earlier experiments to 
larger projectile sizes here. The experiments show that the lim- 
iting filling-factor is independent of the projectile mass when it 
is increased by a factor of 100. Using the same dust material, 
Kothe et al. (2010) found much larger filling factors in multiple- 



collision experiments. ,Kothe et al. (20101 found a filling factor 
of 0.4 at a velocity of 6 ms"' and extrapolated their results to 
values of even 0.6 at 10 ms"'. 

As discussed before, the filUng factor of a growing aggre- 
gate is very important, and the results of |Kothe et al.| ( 20T0| and 
Teiser et al. (2011a) may indicate quite different evolution dur- 



ing further growth. It therefore remains unclear how we can re- 
solve this significant discrepancy, which is one reason why we 
extended the collisions to larger projectile sizes, close to the size 
used by Kothe et"aL] p010). After we had used the same pro- 
jectile masses, only one difference remained in the experiments. 
The targets of Teis er et aL] ( |201 la) l and this work were built by 
multiple impacts at random sites of a target surface, while Kothe] 
|et al.| ( [201 0) studied collisions of always the same site with a 
new aggregate. This created a tip of dust that is smaller than the 
next projectile and is repeatedly compressed. Surface structures 
in our experiments are small compared to the projectile size, and 
the local pressure of an impacting projectile is spread over a 
larger surface area. We propose that this difference in setup as 
an explanation of this discrepancy. Both situations might occur 
in protoplanetary disks. If very irregular aggregates with sharp 
extensions (i.e. from earlier catastrophic disruptions) were hit 
by the next projectile, local filling-factors might be as high as 
determined by Kothe et al. ( 2010| l. However, if random impacts 
created rather smooth surfaces, i.e. a growing aggregate moving 
through a cloud of small aggregates, a lower filling factor with a 
limiting value of about (D = 0.3 1 would result. 

Following up on this are the experiments where aggregates 
were launched at higher velocities up to 6 ms"' but onto a highly 
porous dust target formed in low speed collisions. The experi- 
ments show that the existing part of the aggregate is rather in- 
sensitive to these impacts. Only the new layer that builds up 
consists of dust with a filling factor characteristic of the im- 
pact speed, i.e. that is compressed to a maximum filling-factor of 
0.31. The original aggregate can keep its original high porosity. 
In this way, as impact speed increases with growing target size 
in protoplanetary disks, dust aggregates may consist of a porous 
core and a compact shell. This influences the particle-gas cou- 
pling time, which determines the collision velocities of subse- 
quent collisions, and will likely be important for collisional out- 
comes at larger impact speeds at later phases where collisional 
effects reach the core. Geretshauser et al^ ( |2011b) showed that 
the mechanical properties of these shell-core-aggregates change 
with the ratio of core to shell mass, but is almost dominated by 
the properties of the outer shell. 



& Urbassek (2012l are consistent to within an order of magni- 
tude, considering the different properties of the underlying dust 



Krause et al. ( 201 l| l determined the thermal conductivity of 



dust aggregates in experiments. They varied the filling factor 



8 



T. Meisner et al.: Experiments on centimeter-sized dust aggregates 



and found that more porous dust-aggregates have lower ther- 
mal conductivities. In addition, they found mm-size particles 
drift through a protoplanetary disk at low speed mj s) but 
that this speed still varies significantly across this size range. 
Our experiments show that the filling factor of growing aggre- 
gates is still below the limiting filling-factor of 0.31. As it in- 
creases with velocity and as larger particles are faster, the fill- 
ing factor increases with the size of the aggregates. This, means 
in turn, that larger aggregates have higher thermal conductivity. 
[Loesche & Wurm| ( |subm itted') used these results to estimate the 
photophoretic strength of dust-mantled chondrules, and showed 
that this evolution of filling factors allows a sorting of the parti- 
cles according to size as found in primitive meteorites (KuebleTj 



et al.||1999i. Our experiments show that this is a reaUstic esti- 



mate. 

The experiments and their applications to planet formation 
show that the filling factor of dust aggregates is one of the 
most important parameters determining the physical evolution 
of small bodies in protoplanetary disks. We found that smaller 
(cm-size) aggregates are unlikely to be compact as they could 
be, but eventually reach a limiting filling-factor of 0.31 and have 
a core-shell structure. Since the more compact aggregates are 
more robust in colhsions, they can grow further by feeding on 



the less co mpact aggregates. In e arUer experiments, Wurm et al. 
( 2005b I and |Teiser et al. ( 201 Ib| l demonstrated that larger aggre- 



gates with a filling factor of 0.3 1 can grow even by the impact of 
high speed projectiles. Windmark et al. (2012) showed that this 
process allows the collisional growth of planetesimals. These 
earlier collision experiments used an arbitrary filling-factor of 
0.31 for the manually prepared targets and projectiles. We have 
found and confirmed that this choice was well made as the lim- 
iting filling-factor found in our experiments of self-consistent 
impact growth is exactly the same filling-factor as used in those 
earlier experiments. We conclude that this is strong support for 
a collisional growth scenario of planetesimals. 
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